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ABSTRACT. We study the evolution of the halo-halo correlation function 
and bias in a ACDM model using very high-resolution TV-body simulations 
with dynamical range of ~ 32, 000 (force resolution of ~ 2h~ x kpc and particle 
mass of fa 1O 9 /i _1 M ). The high force and mass resolution allows dark matter 
(DM) halos to survive in the tidal fields in high-density regions and thus 
prevents the ambiguities related with the "overmerging problem." Numbers 
of galaxy-size halos in cluster-like objects in our simulation are similar to 
the numbers of galaxies observed in real clusters. This allows us to estimate 
for the first time the evolution of the correlation function and bias at small 
(down to ~ 100/i _1 kpc) scales. We compare particle distribution, dark matter 
correlation function, density profiles, and halo mass function produced with 
our TV-body code and corresponding results of the AP 3 M simulations. 
We find that at all epochs the 2-point correlation function of galaxy-size halos 
ghh is well approximated by a power-law with slope as 1.6 — 1.8. The difference 
between the shape of £hh an d the shape of the correlation function of matter 
results in the scale- dependent bias at scales < 7h~ 1 Mpc, which we find to be 
a generic prediction of the hierarchical models, independent of the epoch and 
of the model details. We find that our results agree well with existing cluster- 
ing data at different redshifts, indicating the general success of the picture of 
structure formation in which galaxies form inside the host DM halos. Particu- 
larly, we find an excellent agreement in both slope and the amplitude between 
£,hh{ z = 0) m °ur simulation and the galaxy correlation function measured 
using the APM galaxy survey. At high redshifts, the observed clustering of 
the Lyman-break galaxies is also reasonably well reproduced by the models. 



1 Introduction 

It is widely believed that the dis- 
tribution of galaxies is different from 
the overall distribution of dark mat- 
ter (DM). This difference, the bias, is 
crucial for comparisons between ob- 
servations and predictions of cosmo- 
logical models. Observations provide 
information about the distribution of 
objects such as galaxies and galaxy 
clusters. The models, however, most 
readily predict the evolution of the 
dark matter distribution, which can- 
not be observed directly. The mod- 
els, therefore, should be able to pre- 
dict the distribution of objects or, 
conversely, predict how this distribu- 
tion is different from that of the dark 
matter (i.e., the bias). In the seminal 
paper, Kaiser (1984) introduced the 
notion of bias to explain the strong 
clustering of galaxy clusters. Dur- 



ing the subsequent years, many an- 
alytical and numerical studies have 
been done on the subject of the clus- 
ter and galactic bias (see, e.g., refer- 
ences in Colin ct al. 1998). However, 
only recently the dynamic range of 
the pure TV-body simulations has be- 
come sufficiently high to allow to 
study unambigously the small-scale 
bias of DM halos. In this contri- 
bution, we present the results on 
the small-scale (f» 0.1 - 7h~ l Mpc) 
halo-halo clustering and bias of ha- 
los in one of the currently most suc- 
cessful cosmological scenarios: flat 
CDM model with non-zero cosmolog- 
ical constant (ACDM). The parame- 
ters of the model are as follows: f2o = 
1 - fi A = 0.3; h = 0.7; ct 8 = 1-0; 
and age of ss 13.4 Gyrs. These re- 
sults are a part of a larger study of 
the bias evolution in different cosmo- 
logical models described in Colin et 
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Figure 1. The comparison of the cor- 
relation functions of the dark mat- 
ter £dm in the ACDM model es- 
timated by different authors with 
different numerical resolutions and 
codes. The solid curve shows £d m 
in our ACDM run, simulated us- 
ing the ART code. The dashed curve 
shows estimated by Jenkins et 

al. (1998, ApJ, 499, 20) using the 
AP 3 M code. The crosses show 
estimated by Klypin et al. (1996, 
ApJ, 466, 13) using the PM code. 
The vertical lines indicate formal 
force resolution for each code (the 
line for the ART code at l.Sh -1 kpc 
is off the plot). 

al. (1998) and we refer the reader to 
this paper for further details. 

2 Numerical Simulation 

One of the main goals of our 
study was to compute the corre- 
lation function and the bias ac- 
counting for all DM halos, includ- 
ing those inside groups and clusters. 
To assure that halos do survive in 
clusters the force resolution should 
be ~ 1 — 3/i~ 1 kpc (Moore et al. 
1996; Klypin ct al. 1998). Further- 
more, if we aim to study galaxy- 
size halos, the mass resolution should 
be < 1O 9 /i _1 M0 to resolve galaxy- 
size halos (M > 10 11 /i- 1 M Q ) with at 
least f» 100 particles. The compro- 
mise between these considerations 
and the computational expense de- 
termined the force and mass res- 



olution of our simulation carried 
out with the Adaptive Refinement 
Tree (ART) N-body code (Kravtsov, 
Klypin & Khokhlov 1997): particle 
mass of Ri 1O 9 /i -1 M0 and peak spa- 
tial resolution of r//2 Ri 2h~ 1 kpc. 
The dynamical range Li, ox /r] of the 
simulation implied by the force res- 
olution is £3 16,000 (32,000 formal). 
The ART code integrates the equa- 
tions of motion in comoving coordi- 
nates. However, its refinement strat- 
egy is designed to effectively pre- 
serve the initial physical resolution 
of the simulation. In order to pre- 
vent degradation of force resolution 
in physical coordinates, the dynamic 
range between the start and the end 
(z = 0) of the simulation should 
increase by (1 + Zj): i.e., for our 
simulations reach 512 X (1 + Zj) = 
15, 872. This is accomplished with 
the prompt successive refinements 
in high- density regions during the 
simulations. The peak resolution is 
reached by creating a refinement hi- 
erarchy with six levels of refinement. 
The spatial refinement is accompa- 
nied by the similar refinement of in- 
tegration time step. 

We use Bound Density Maxima 
(BDM, Klypin et al. 1998) algo- 
rithm to identify galaxy-size halos 
that are isolated or belong to a 
larger system. The main idea of the 
BDM algorithm is to find positions 
of local maxima in the density field 
smoothed at a certain scale and 
to apply physically motivated cri- 
teria to test whether the identified 
site corresponds to a gravitationally 
bound halo. The detailed description 
of the halo finder and simulation pa- 
rameters can be found in Colin et al. 
Here we present results of the tests of 
code performance based on compar- 
isons with AP 3 M simulations. In Fig- 
ure 1 we compare the 2-point corre- 
lation function of the dark matter in 
our simulation with similar estimates 
made with the PM and AP 3 M codes 
(cf. fig. caption). We have made ap- 
propriate corrections to the ampli- 
tude of the correlation function to 
account for slightly different normal- 
izations of the simulations. Figure 
1 shows that there is a very good 
agreement between all estimates at 
scales ss (0.2 -2)h- 1 Mpc. The ART 
and AP 3 M estimates agree to better 
than 10% at scales 0.03 -Ih' 1 Mpc! 
On larger scales the ART correla- 
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Figure 2. Comparison of z = outputs of ART and AP^M CDM simulations 
started from the same initial conditions. Two bottom panels compare particle 
distributions in a slice through the simulation box. The upper panels compare 
halo density profiles (left) and halo mass functions (right). The dashed (solid) 
lines show the ART (AP^M) results. 



tion function has a smaller ampli- 
tude than that of the AP 3 M due to 
a factor of 2.35 smaller box size of 
the former. This result is in agree- 
ment with previous studies which 
concluded that the correlation func- 
tion is underestimated at scales > 0.1 
of the simulation box size. Never- 
theless, the agreement is striking 
at smaller scales, given all the dif- 
ferences (including cosmic variance) 
between the estimates. We conclude, 
therefore, that the correlation func- 
tions discussed in the next section 
are robust at scales < 7h~ 1 Mpc. 

The second test used smaller sim- 
ulations of the CDM model to test 
the distribution and properties of ha- 
los produced by the ART and AP 3 M 
codes. This test is a part of a com- 
parison study carried out in collab- 
oration with Alexander Kncbc and 
Stefan Gottlober (AIP, Potsdam). 
The ART simulations had a (for- 
mal) dynamic range of 4096, while 



the AP 3 M simulation had a dynamic 
range of 6400. Figure 2 shows com- 
parison of the particle distributions, 
density profiles of halos, and halo 
mass functions in the two simula- 
tions. Although the small-scale de- 
tails of the particle distribution are 
slightly different (due to different 
time integration schemes used and 
somewhat different spatial resolu- 
tions) the overall distribution is very 
similar in the two runs. Almost iden- 
tical results are also observed for the 
density profiles and halo mass func- 
tion. 



3 Results 

Our analysis shows that at all 
epochs correlation functions of ha- 
los and dark matter arc different in 
both shape and amplitude at scales 
< 7/i — 1 Mpc. The difference in am- 
plitude means that halos are biased 
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with respect to dark matter, while 
the difference in shape implies that 
the bias is scale- dependent. Figure 3 
shows evolution of the bias with red- 
shift at different scales and for halos 
with different minimum circular ve- 
locities^]. We use a conventional defi- 
nition of the bias as a square root of 
the ratio of the halo and dark mat- 
ter correlation functions. The figure 
shows that bias undergoes a signif- 
icant evolution with redshift: it is 
high, b ~ 2 — 4 at high redshifts 
z ~ 3 — 5, and then rapidly decreases 
at lower redshifts to b ~ 0.5 — 1 at 
z = 0. The high value of bias at early 
epochs is in reasonably good agree- 
ment with the recent observations of 
strong galaxy clustering at high red- 
shifts (Steidel et al. 1998). The Fig. 

3 also shows that at high redshifts, 
bias depends on the mass of halos: 
the massive halos are clustered more 
strongly than less massive ones. The 
scale-dependency of the bias is illus- 
trated in the bottom panel of the 
Fig. 3, in which we show evolution of 
the bias at three different scales for 
the halos with the circular velocity > 
lOOkm/s. Although qualitatively the 
evolution is similar, the rate of evolu- 
tion and the value of bias are clearly 
dependent on the scale. In this panel 
we also show analytical prediction of 
given by Moscardini et al. (1998) for 
halos of mass > 1O 11 /i~ 1 M . There 
is a good qualitative agreement be- 
tween analytical predictions and our 
numerical results. 

We find that at all epochs the cor- 
relation function of DM halos is well 
described by a power law £(r) = 
(r/ro) 7 , where ro ~ 3 — 5h~ 1 Mpc 
and 7 rj —1.7. This is in good agree- 
ment with the shape and amplitude 
of the galaxy correlation function at 
both low and high redshifts. Figure 

4 shows the z = correlation func- 
tions of halos in our simulation and 
the correlation function of galaxies 
from the APM galaxy survey. The 
figure shows that there is an excel- 
lent agreement between the two at 
all scales probed by our simulation. 
The dotted line in the bottom panel 
of Fig. 4 shows the correlation func- 
tion of the dark matter. Note the dif- 
ference in shape between correlation 

1 The circular velocity is a mass in- 
dicator: the more massive is a halo, 
the higher is its circular velocity. 
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Figure 3. Top panel: The evolu- 
tion of bias at comoving scale of 
0.54h -1 Mpc for halos with differ- 
ent lower limit on the maximum cir- 
cular velocity in the ACDM simula- 
tion. The bottom panel shows depen- 
dence of the bias on (comovingj scale 
for halos with maximum circular ve- 
locity > lOOkm/s. 

functions of halos and dark matter. 
This difference implies that the bias 
is scale-dependent; the bias as a func- 
tion of scale is shown in the upper 
panel of Fig. 4. Note also that halo 
correlation function has a lower am- 
plitude than that of dark matter at 
scales < 5?i — 1 Mpc. This means that 
halos at these scales are anti-biased 
with respect to the dark matter. 

Overall, the comparisons dis- 
cussed above and more detailed anal- 
ysis of Colin et al. (1998) indicate 
that there is good agreement be- 
tween our results and the clustering 
data at both low and high redshifts. 
This implies that hierarchical mod- 
els in which observed galaxies form 
in the host DM halos naturally ex- 
plain the observed galaxy clustering 
at different epochs, including excel- 
lent agreement with the accurately 
measured 2 = correlation function. 
On the other hand, the generic form 
of the bias evolution observed in the 
numerical simulations at high red- 
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Figure 4. Bottom panel: Comparison of the halo correlation functions (solid, 
dot- dashed, and dashed curves) in the simulation with the correlation function 
of the APM galaxies (circles; Baugh 1996, MNRAS, 280, 267). The dotted 
curve shows the dark matter correlation function. Top panel: Dependence of 
bias on scale and maximum circular velocity. The curve labeling is the same 
as in the bottom panel, except that the dotted line now represents the bias of 
halos with V m ax > lOOkm/s. 



shifts agrees well with the prediction 
of the analytical models based on the 
extended Press-Schechter formalism. 
This implies that we understand the 
nature of the bias and the processes 
that drive its evolution at high z. 
It is also reassuring that our results 
are in good agreement with results 
of the studies which have used semi- 
analytical modelling to follow evolu- 
tion of galaxies in DM halos (e.g., 
Kauffmann et al. 1998; Benson et al. 
1998), as well as with the results of 
direct numerical hydro simulations 
(Jenkins et al. 1998; this meeting). 
The success of the current theoretical 
models in interpreting the cluster- 
ing data forms a solid foundation for 
further sophistication of the numeri- 
cal models by including the processes 
important for galaxy formation (such 
as dynamics of baryons, cooling, star 
formation, and stellar feedback). We 
can expect, therefore, that the prob- 
lem of galaxy biasing will be tack- 
led in the foreseeable future, which 
would open for us great possibilities 
for the analysis of the ever increas- 



ing amount of observations of galaxy 
clustering at low and high redshifts. 
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